Spectrum of cosmic rays with energy above 10 eV 
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There are some discrepancies in the results on energy spectrum from Yakutsk, AGASA, and HiRes experiments. 
In this work differential energy spectrum of primary cosmic rays based on the Yakutsk EAS Array data is presented. 
For the largest events values of S^oo and axes coordinates have been obtained using revised lateral distribution 
function. Simulation of converter's response at large distances showed no considerable underestimation of particle 
density. Complex shape of spectrum in region of Eo > 10 17 eV is confirmed. After adjustment of parameters and 
additional exposition at the Yakutsk array there are three events with energy Eo > 10 20 eV. 



1. INTRODUCTION 

Research into the cosmic rays above 10 17 eV 
spectrum shape and into intensity in steep re- 
gion near 10 20 eV predicted by Greisen pQ, Zat- 
sepin and Kuzmin [5] are of great importance 
for their sources detection. Results obtained in 
various experiments |3I4I5| differ from each other 
by factor 2 and more in absolute intensity, but 
their shapes are similar. Variations of intensity 
cannot be described with simple power law. At 
energies greater than the GZK-cutoff results are 
inconsistent. At HiRes there is only one event 
with Eq > 10 20 eV and the spectrum is cut-off. 
AGASA has registered 11 such events (6 < 45°), 
this could be an evidence for absence of a cutoff. 

At the Yakutsk array after recent analysis have 
been carried out |3] there is only one event with 
energy estimated to be greater than 10 20 eV. To 
explain this contradiction with AGASA A. Wat- 
son assumed that at the Yakutsk array such show- 
ers are skipped due to inadequate short integra- 
tion time for large distances from the axis. In 
this work we have studied affect of particle ar- 
rival time distribution at different distances on 
estimated density for Yakutsk and AGASA. Sim- 
ulation showed that distortion of estimated den- 
sity at the Yakutsk array cannot result in such 
miscalculation. More substantively is the lateral 
distribution function (LDF) used for axis deter- 



mination. We provided axes coordinates determi- 
nation with adjusted LDF for the largest events, 
which in average led to increase of S^oo- Be- 
sides, February 18 2004 a new event with energy 
10 20 eV was detected at the Yakutsk array. As 
a result of these factors there are three events of 
Eo > 10 20 eV registered at the Yakutsk array. 

2. DENSITY MEASUREMENT AT 
LARGE DISTANCES FROM THE 
AXIS 

A shower disk at lagre distance from the axis 
(R > 1500 m) is quite thick. At Yakutsk ar- 
ray and at AGASA a nearly similar logarithmic 
RC-converters of the signal from photomultiplier 
tube (PMT) to digital code with r ~ 10 mcsec 
are used. At the Yakutsk array for an event to be 
treated, a coincidence of signals from both detec- 
tors within 2 mcsec is required. Herewith input of 
converters is closed in 2 mcsec after coincidence. 
In the case when shower front is wide, this may 
result in underestimation of the density. Besides, 
there is a possibility, that station doesn't operate 
due to large difference in particles arrival times on 
different detectors and so some showers might be 
skipped. These circumstances have been pointed 
out by Watson in his report 6. . At AGASA input 
is permanently open and in the case of wide sig- 
nal this may lead to density overestimation due 
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to features of converter. 

To examine the influence of the effects men- 
tioned above, we have provided simulation of de- 
tector's response for distances R = 1050, 1500 
and 2000 m, based on the particle distribution 
approximation obtained at AGASA [7]. A coeffi- 
cient Kr was considered — a ratio between den- 
sity estimated with RC-convector and the one 
set with program. For the Yakutsk array, in 
the case of detectors with large area, register- 
ing large particle densities we got following val- 
ues: ifioso = 1-05, K 150 o = 0.994, K 20 oo = 
0.76. Same points for AGASA: K W5Q = 1.065, 
#1500 = 1-11, #2000 = 1-2. At 2000 m dis- 
tance for Yakutsk there is 25% underestimation, 
for AGASA — 20% overestimation. 

For the real experiment for the shower with 
E = 10 20 eV at R = 2000 m about 2 parti- 
cles per detector is expected. Simulation indi- 
cated, that in this case underestimation is much 
less than #2000 = 0.92. It is connected with the 
fact that conversion starts only after the first par- 
ticle hit and at low density the effective thickness 
of the shower front decreases. Probability of that 
the station doesn't operate due to gap between 
operating of two separate detectors is more than 
2 mcsec is 8.5% and it is lower by factor 3 than 
those due to Poisson fluctuations at this density. 
But considering this circumstance together with 
steepness of the LDF at large distances we re- 
stricted used effective area beyond the bound of 
array to obtain the intensity of the largest show- 
ers. 

Simulation showed no significant underestima- 
tion of particle density for distances up to 2000 m 
for density estimating system at the Yakutsk ar- 
ray. In the case of AGASA, when input of RC- 
convector is constantly open, besides wide distri- 
bution, there is an afterpulse contribution to den- 
sity overestimation from delayed particles (proba- 
bly neutrons) together with casual additives from 
background muons. One can conclude from the 
data in |7] that lagging neutrons can overstate 
the density by 1.37 already at 500 m and further. 
Background muons may cause distortions in wide 
range of axis distances. If one such particle hits 
within last 10 mcsec of RC-circuit discharge, then 
resulting density can be overestimated by factor 2 



and more independently of the real density. The 
effects mentioned above are excluded thanks to 
closing of converter's input in 2 mcsec. 

3. ENERGY SPECTRUM SUMMARY 

Events selection was provided as described 
in PJ. Showers with 9 < 60° were used. For 
determination of the intensity for showers with 
Eq > 4 • 10 19 eV, an extended area together effi- 
cient zone outside the array was used. 

At standard procedure of axis determination, 
the Greisen-Linsley approximation of LDF with 
parameters obtained before at the Yakutsk ar- 
ray is used [Hj- It was shown in that for 
showers with energy greater than 10 19 eV this 
LDF badly corresponded with experimental data 
at R > 1000 m from the axis. A modified approx- 
imation was proposed: 
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For showers with Eq > 2 ■ 10 19 eV parameter b 
does not depend on energy but depends on 9. In 
this work we provided axes coordinates determi- 
nation with this adjusted LDF. As a result, in 
average estimated S600 values increased: 10% — 
for showers with axes lying within array area and 
20% on border. 

For energy estimation we used adjusted formu- 
las of 5300 and S^oo E3- F° r upright showers 
(6> = 0°, depth X Q = 1020 g ■ cm" 2 ), the follow- 
ing relationships were obtained using calorimetric 
method: 
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E = (4.6 ± 1.2) • 10 17 x 
x [W0°)]°- 98±0 - 02 eV. 
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Zenith-angular dependency was accepted in 
form of: 
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where A c = 200 g/cm 2 - - the attenuation 
length for the soft component (electrons), A m = 
1000 g/cm 2 — the attenuation length for the hard 
component (associated with muons), (3 — contri- 
bution from the hard component to full response 
5 3 oo(0°) or S 60 o(0°) at 1020 g ■ cm" 2 depth. 

The following dependencies of parameters p on 
5300 and Seoo have been carried out from experi- 
mental data: 
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Figure 1. Differential energy spectra accord- 
ing to the Yakutsk EAS array, AGASA g] and 
HiRes 0. 



are presented. Results obtained in different ex- 
periments correspond quite well in shape but dif- 
fer in intensity. The data from the Yakutsk ar- 
ray near 10 19 eV are higher by factor ~ 2.5 than 
HiRes data and - 30% than AGASA's. This 
rather is connected with the difference in esti- 
mation of the showers energy. For estimation of 
energy in upright showers at AGASA the depen- 
dency between 5*600 and Eq based upon the model 
calculations [7| is used. Our formulas J2J and J3J 
for upright showers estimate the energy (30-40)% 
higher at E « 5 • 10 17 eV and (15-20)% higher at 
Eq ~ 10 19 eV than it follows from similar models. 

In the region of ultra-high energies the re- 
sults from Yakutsk and AGASA approach. There 
are 4 events registered in Yakutsk with adjusted 
estimated energy exceeding 10 19 9 eV (see Ta.- 
ble^l. Relative errors in energy estimation re- 
sulting from uncertainty of the parameters in for- 
mulas and errors in determination of axis coordi- 
nates and zenith angle in these individual show- 
ers amount from 32% to 46%. If the energy is 
reduced by one standard error then it slightly ex- 
ceeds the 10 20 eV threshold only in one event. 
Therefore the relic cutoff of the spectrum cannot 
be rejected based on Yakutsk EAS data. 

Similar experimental errors are observed at 
AGASA. According to ^H] their averaged value 
is about 20%. Taking into account this circum- 
stance a conclusion was made in ,11 that yet 
there are too few events recorded to approve the 
spectrum cutoff absence. Besides, estimations of 
the energy at AGASA depend on model conclu- 
sions. The affects of observed densities mentioned 
above are also not considered yet. 

The HiRes results are consistent with the 
GZK-cutoff of the spectrum, the AGASA and 
Yakutsk data are inconsistent. But because of 
small statistics and errors in energy estimation 
while it is impossible to final conclude about this 
problem. To solve this and investigate of affinities 
of the particles with energies above GZK-cutoff, 
data with high statistics and good accuracy in the 
energy estimation a necessary available. 



On Fig^ the differential energy spectra ob- 
tained at the Yakutsk array, AGASA and HiRes 
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Table 1 



The most energ 


etic events detected with the Yakutsk array (sorted by 


energy) 




N Date 


Time, UT 


6° 


log So 


SE (%) 


b° 


r 


1 02.18.04 


22:20:38 


47.7 


20.16 


42 


16.3 


140.2 


2 05.07.89 


22:03:00 


58.7 


20.14 


46 


2.7 


161.6 


3 12.21.77 


18:45:00 


46.0 


20.01 


40 


50.0 


220.6 


4 02.15.78 


03:35:00 


9.6 


19.99 


32 


15.5 


102.0 
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